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ABSTRACT 

Based on the results of applying the extended ADC emission model for low mass X-ray binaries to three Z-track sources: 
, GX 340-1-0, GX5-1 and CygX-2, we propose an explanation of the Cygnus X-2 like Z-track sources. The Normal Branch is 

■ dominated by the increasing radiation pressure of the neutron star caused by a mass accretion rate that increases between the 
' soft apex and the hard apex. The radiation pressure continues to increase on the Horizontal Branch becoming several times 
, super-Eddington. We suggest that this disrupts the inner accretion disk and that part of the accretion flow is diverted vertically 

■ forming jets which are detected by their radio emission on this part of the Z-track. We thus propose that high radiation pressure 
is the necessary condition for the launching of jets. On the Flaring Branch there is a large increase in the neutron star blackbody 

■ luminosity at constant mass accretion rate indicating an additional energy source on the neutron star. We find that there is 
' good agreement between the mass accretion rate per unit emitting area of the neutron star m at the onset of flaring and the 

theoretical critical value at which burning becomes unstable. We thus propose that flaring in the Cygnus X-2 like sources consists 
of unstable nuclear burning. Correlation of measurements of kilohertz QPO frequencies in all three sources with spectral fitting 
results leads to the proposal that the upper kHz QPO is an oscillation always taking place at the inner accretion disk edge, the 
, radius of which increases due to disruption of the disk by the high radiation pressure of the neutron star. 

Key words. Accretion: accretion disks - acceleration of particles - binaries: close - stars: neutron - X-rays: binaries - X-rays: 
. individual (CygX-2, GX 340-^0, GX5-1) 



1. Introduction 

The Z-track sources are the brightest group of Galactic 
low mass X-ray binaries (LMXB) containing a neutron 
star persistently emitting at the Eddington luminosity or 
several times this. The sources trace out a Z-shape in 
hardness-intensity or hardness-softness (Hasinger & van 
der Klis 1989) having three branches labelled the hori- 
zontal branch (HB), the normal branch (NB) and flar- 
ing branch (FB). Strong physical changes are clearly tak- 
ing place at the inner disk and neutron star, but a gen- 
erally accepted explanation of the Z-track phenomenon 
has not existed. Hasinger & van der Klis (1989) presented 
colour-colour diagrams for all six original Z-track sources 
using Exosat data which showed rather different shaped Z- 
tracks between two groups: Cyg X-2, GX 340-f-O and GX 5- 
1 (the Cygnus X-2 like sources) and ScoX-1, GX 17-1-2 and 
GX 349-1-2 (the ScoX-1 like sources). It was proposed that 
these constituted two sub-groups of Z-track sources, for 
example by Kuulkers et al. (1994), and this became gen- 
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erally accepted (Kuulkers & van der Klis 1995; Smale et al. 
2000, 2003). Many observations have confirmed the rather 
different shaped Z-tracks, the Cygnus X-2 like sources hav- 
ing a long HB, while the ScoX-1 like sources have a short 
HB but a long FB. 

LMXB have been classified according to the various 
phenomena observed, into the Atoll, Z-track, dipping and 
ADC sources, and this has been shown to relate to in- 
clination angle so that X-ray dipping is seen at inclina- 
tions greater than ^65°. The Z-track sources display char- 
acteristic patterns in hardness-intensity but, in general, 
no orbital related behaviour indicating smaller inclina- 
tion and the Atoll sources exhibit rather different patterns 
of hardness-intensity. However, a more basic classification 
would be based on luminosity in which the Z-track sources 
have luminosities greater than ~10^^ erg s~^ and the Atoll 
sources have luminosities between ^10"^^ and ^10"^^ erg 
s^^, although not all lower luminosity LMXB have been 
identified as Atoll sources via their hardness- intensity pat- 
terns. These patterns are also not understood, nor is the 
relation between the Z-track and Atoll classes, so that 
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understanding Z-track and Atoll sources remains a funda- 
mental problem. 

Thus we need to understand the detailed nature of the 
inner disk and accretion in the super-Eddington LMXB, 
the Z-track sources, but these are also important since 
they are detected as radio emitters, although essentially 
in one part of the Z-track only: the horizontal branch. Not 
only is radio detected, but striking results from the VLA 
show the release of massive radio condensations from the 
source ScoX-1 (Fomalont et al. 20001) moving outwards 
in two opposite direction at velocity v/c ~ 0.45. Because 
radio is detected essentially in one branch only, the sources 
offer the possibility of determining conditions pertaining 
inside the sources via X-ray observations, i.e. the condi- 
tions at the inner disk, when jets are present, which dis- 
tinguish the horizontal branch from the other branches, 
and so finding the conditions necessary for jet formation. 

Possible ways of understanding the Z-track sources are 
by timing studies, spectral studies or theoretical studies. 
Extensive timing studies have been made to investigate 
the quasi periodic oscillations (QPO) found in the Z-track 
and other sources (e.g. van der Klis et al. 1987). This has 
revealed a wealth of complex behaviour, and in particular, 
study of the kilohertz QPO has revealed a pattern of vari- 
ation of the twin kHz peaks on the horizontal branch of 
Z-track sources, and various models have been proposed 
to explain the QPO (Stella & Vietri 1998; Abramowicz &: 
Kluzniak 2001; Miller et al. 1998). However, it is difficult 
to determine what changes take place at the inner disk 
around the Z-track from data on the kHz QPO. 

Previous spectral fitting tended to apply a particular 
LMXB emission model: the Eastern model of Mitsuda et 
al. (1989) in which the emission is assumed to consist of 
disk blackbody emission plus Comptonized emission from 
the neutron star or inner disk. Application of this model 
may have given acceptable fits to the spectra (Done et al. 
2002, Agrawal & Sreekumar 2003; di Salvo et al. 2002), 
however, interpretation of the results, i.e. of the variation 
of the various spectral parameters around the Z-track was 
not easy. Our previous work over an extended period on 
the dipping class of LMXB provided evidence for a differ- 
ent emission model and we proposed the "extended ADC 
(accretion disk corona)" model (Church & Baluciriska- 
Church 2004, 1995) in which the emission components 
consist of blackbody emission from the neutron star plus 
Comptonized emission from an extended ADC. In this 
model, the ADC covers a substantial fraction of the in- 
ner disk and the thermal emission of the disk provides the 
seed photons for Comptonization (Church & Balucinska- 
Church 2004). The model was developed from longterm 
studies of the dipping LMXB sources. In all LMXB, the 
observed shape of the spectrum and spectral fitting show 
that Comptonized emission dominates the spectrum, i.e. 
comprises more than 90% of the total luminosity, as shown 
for example in the ASCA survey of LMXB of Church & 
Balucinska-Church (2001). In less luminous sources, the 
percentage may be higher, while in the brightest sources, 
the Z-track sources, the contribution is ~90% at the soft 



apex. In the dipping sources, it is generally accepted that 
the reductions in X-ray intensity at the orbital period are 
due to absorption in the bulge in the outer accretion disk 
(White & Swank 1982; Walter et al. 1982). Spectral fitting 
of non-dip and dip spectra clearly shows that the domi- 
nant Comptonized spectral component is removed only 
slowly in dipping and so the 90% contributor to the to- 
tal luminosity must be extended (e.g. Church et al. 1997). 
Measurement of the dip ingress times provides the sizes 
of the extended emission. These measurements indicated 
the radial extent to be very large: between 20 000 and 
700 000 km depending on source luminosity (Church & 
Balucinska-Church 2004). Moreover, this extended region 
cannot reach very high above the disk otherwise it would 
not be possible to fully cover the emitter by absorber and 
100% deep dipping, often observed, would not be possible. 
Thus the conclusion cannot be avoided that 90% of the 
total luminosity originates in this extended region which, 
given the measured geometry: large in radial extent but 
small in height, can be identified with a hot, thin accretion 
disk corona above the disk, typically extending across 15% 
of the disk radius. Recently independent support for the 
model came from the Chandra grating results of Schulz et 
al. (2009) on CygX-2. Precise grating measurements re- 
vealed a wealth of emission lines of highly excited states 
such as the H-like ions of Ne, Mg, Si, S and Fe. The width 
of these lines indicated Doppler shifts due to orbital mo- 
tion in the accretion disk corona at radial positions be- 
tween 18 000 and 110 000 km in good agreement with the 
overall ADC size from dip ingress timing. Further support 
for line formation at large radial distances in an extended 
ADC came from recent resiilts of Chandra observations 
of 4U 1624-490 (Xiang et al. 2009). The evidence for the 
extended corona is thus now overwhelming. 

Theoretical modelling for the Z-track soiirces was car- 
ried out by Psaltis et al. (1995) assuming a magnetosphere 
of the neutron star and the changing properties and geom- 
etry of this as the mass accretion rate changcxl. Howc!V(?r, 
the model assumed that the Comptonized emission ob- 
served in the spectra of all LMXB originated in a small 
central region close to the neutron star, not in an extended 
region. 

It should also be stated here that although the na- 
ture of the Z-track has not been clear, there has been a 
standard assumption made that the mass accretion rate 
increases monotonically around the Z-track in the sense 
HB - NB - FB, based originally on results from multi 
wavelength observations of CygnusX-2 (Hasinger et al. 
1990). This assumption has been widely but not univer- 
sally adopted and it will appear in Sects. 3 and 4 that our 
results do not support this assumption. 

More recently, we have taken the approach of applying 
the extended ADC model to the Z-track sources GX 340-1-0 
(Church et al. 2006) and GX5-1 (Jackson et al. 2009), on 
the basis that the evidence favours the Extended ADC 
model, to test the hypothesis that this can fit the spec- 
tra, and provide a physically reasonable and consistent 
explanation of the Z-track. Good fits to the spectra at 
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Fig. 1. Background-subtracted and deadtime-corrected PCA light curve of CygnusX-2 with 64 s binning. 



all positions on the Z-track were obtained, and the spec- 
tral fitting clearly suggested an explanation of the physical 
changes taking place on the Z-track. At the Soft Apex: the 
junction of the HB and FB, the X-ray intensity is mini- 
mum. Ascending the normal branch to the Hard Apex, 
the intensity increases substantially and the neutron star 
blackbody temperature kT increases from ~1 to ~2 keV. 
The luminosity of the ADC Comptonized emission (Ladc) 
which is 90% of the total 1-30 keV luminosity at the soft 
apex increases by a factor of two or more strongly suggest- 
ing that M was increasing. The increase of kT is consistent 
with this and has a major consequence that the radiation 
pressure of the neutron star (~ T*) increases by an or- 
der of magnitude. The strength of the radiation pressure 
close to the neutron star which will determine its effects 
on the inner accretion disk is conveniently expressed in 
terms of the ratio of the flux emitted per unit area to its 
Eddington value ///Edd where f^dd = L-Edd/^T^R^ and 
R is the radius of the neutron star. This ratio is a bet- 
ter measure of radiation pressure effects close to the neu- 
tron star than L/L^dd which averages the flux over the 
surface, especially if only part of the surface is emitting. 
In GX 340+0, ///Edd was found to increase from 0.2 at 
the soft apex to 1.5 at the hard apex rising to almost 3 
times super-Eddington on the HB. In GX5-1 the ratio 
similarly increased to become greater than unity, so that 
in both sources, the radiation became super-Eddington at 
the hard apex and HB. It is these parts of the Z-track on 
which jets are observed and we proposed that the high ra- 
diation pressure disrupts the inner disk and drives matter 
into the vertical direction so launching the jets (Church 
et al. 2006; Jackson et al. 2009). 

In both sources, on the flaring branch Labc was con- 
stant suggesting M was constant, while the neutron star 
blackbody luminosity increased strongly implying that a 
non-accretion powered energy source developed on the 
neutron star. The mass accretion rate per unit area of the 
neutron star to = M/4 tt Rqq, where i?BB is the blackbody 
radius derived from spectral fltting, was found at the start 



of the FB, i.e. the soft apex, to be in good agreement with 
the theoretical critical value msx (Bildstcn 1998). This is 
the critical value of m below which nuclear burning on 
the surface of the neutron star becomes unstable, i.e. for 
TO < 1.3 X 10^ g cm^^ s^^ the burning changes from sta- 
ble to unstable. Thus it was proposed that for a source 
descending the normal branch with falling M, the onset 
of unstable burning causes the strong flaring seen in the 
X-ray lightcurve. The increase in blackbody luminosity is 
simply the nuclear burning power, as seen at essentially 
constant LadC) i-e. constant M. 

In the second source, GX 5-1, the evolution of the kilo- 
hertz QPO was also investigated by Jackson et al. (2009). 
Using the same selections of data used in spectral analy- 
sis, the variation of the kHz QPO along the HB and part 
of the NB was obtained by timing analysis, and viewed 
in the light of the spectral fitting results. It is well-known 
that the upper and lower kHz QPO frequencies 1/2 and 
i^i vary on the HB (e.g. Jonker et al. 2000; Wijnands et 
al. 1998a,b), and that if U2 is assumed to correspond to 
a Keplerian frequency in the disk, the radial position r 
is ~18 km, i.e. some way into the disk. A previously un- 
known correlation was found between the measured fre- 
quencies and the above ratio / //Edd (Jackson et al. 2009) , 
specifically that the r was approximately constant until 
///sdd = 1) above which r increased rapidly. This sup- 
ports the proposal of disk disruption by high radiation 
pressure; moreover it suggests a mechanism for the upper 
frequency kHz QPO that this oscillation is always an os- 
cillation at the inner disk edge, the radius of which was 
modified by different levels of radiation pressure. Thus 
disruption of the inner disk was the factor determining 
the frequency of the higher frequency kHz QPO and its 
variation. 

In this work, we apply the extended ADC model for 
the first time to the Z-track source CygnusX-2, also deter- 
mining the evolution of the kHz QPO around the Z-track 
in Cyg X-2 and in GX 340-1-0 and we present a comparison 
of the three Cygnus X-2 like sources. 
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2. Observations and analysis 

We analysed the Rossi-XTE ohscivation of CygnusX-2 of 
June 30 - July 3, 1997 in which a full Z-track was covtacxl 
in a relatively short time. The observation consisted of 11 
sub-observations spanning 235 ksec. Data from both the 
proportional counter array (PCA: 2-60 kcV) and the high 
energy X-ray timing experiment (HEXTE: 15 - 250 keV) 
were used. PCA data in Standard2 mode with 16 s resolu- 
tion were used for the extraction of light curves and spec- 
tra. Examination of the housekeeping data showed that 
three of the five Xe proportional counter units (PCUs) 
were consistently on, and so data were extracted from 
these (PCUs 0, 1 and 2). Standard screening criteria were 
applied to select data with an offset between the source 
and tc;lcscope pointing direction of less than 0.02° and 
elevation above the Earth's limb greater than 10°. Data 
were extracted from the top layer of the detectors using 
both left and right anodes. Analysis was carried out using 
the standard RXTE software ftools 5.3.1. A lightcurve 
was generated in the band 1.9 - 20.3 keV and the fa- 
cility PCABACKEST used to generate background files for 
each raw data file based on the latest bright-source back- 
ground model for Epoch 3 of the mission during which the 
observations were made. Deadtime correction was carried 
out, and the background-subtracted, deadtime-corrected 
lightcurve is shown in Fig. 1 with 64 s binning. Lightcurves 
were also extracted in throe sub-bands: 1.9 - 4.1, 4.1 - 6.9 
and 6.9 - 20.3 kcV. By dividing the high band by the 
medium band lightcurve, a hardness-intensity plot was 
made as shown in Fig. 2. Individual subsections of the 
lightcurve were identified with their positions on the Z- 
track by making a hardness- intensity plot for each and lo- 
cating its position on the overall Z-track. Thus the three 
Z-track branches are identified as shown in Fig. 1. 

Lightcurves and spectra were also extracted from 
Cluster 1 and Cluster 2 of the HEXTE instrument, so that 
simultaneous spectral fitting of PCA and HEXTE spec- 
tra could be carried out. These were extracted with the 
FTOOL HXTLCURV which also provided background files 
and allowed deadtime correction. 

Timing analysis was carried out as a function of posi- 
tion along the Z-track. The data used for this were single 
bit data of type SB_125us_14_249_ls, i.e. data having 125 
/xs resolution and a single energy bin spanning 5-60 keV. 
Power density spectra were extracted for each selection on 
the Z-track using the Xronos facility powspec in a single 
energy bin, and Leahy normalization was applied. To allow 
high quality fitting of the spectra, they were converted to a 
format allowing use of the Xspec spectral fitting software, 
by modifying the file headers and re-naming of columns. 
FLx2xsP was used to produce a file with the required pha 
extension and to generate a dummy instrument response 
as needed by Xspec consisting of a unitary matrix having 
no effect on the data. The spectra were grouped to a min- 
imum count per bin to allow use of the statistic, the 
grouping varying between 40 and 140 counts per bin as 
found appropriate by trial and error. 



Timing analysis was also carried out for the source 
GX 340+0 (Sect. 5) so that in all three sources timing 
analysis results could be correlated with spectral fitting 
results for the same selections of data along the Z-track. In 
the case of GX 340+0, the data used were event mode data 
of type E_8us_32B_14_ls which were fitted in the range 5 
- 60 keV. The data were regrouped appropriately, again 
having 4 - 140 counts per bin. 

Following presentation of results for CygnusX-2, we 
present a comparison of spectral fitting and of timing 
analysis results for all three CygnusX-2 like sources: 
CygX-2, GX 340+0 and GX5-1 using our previous analy- 
ses of the RXTE observations of GX 340+0 made on 1997 
September lasting 400 ksec and the 1998 November ob- 
servation of GX5-1 spanning 95 ksec (Church et al. 2006; 
Jackson et al. 2009). Results for the 1997 June observation 
of GX 340+0 are also shown. 

3. Results: CygnusX-2 

The PCA lightcurve shown in Fig. 1 shows that in the 
first 25 ksec of the observation, the source was in a period 
of fiaring, moving on the Z-track along the flaring branch, 
and then onto the normal branch. The source stayed on 
the normal branch, until at ~100 ksec it moved to the 
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horizontal branch. It can be seen that at about 200 ksec 
the intensity of the source had drifted upwards in compar- 
ison with the previous parts of the observation. Plotting 
these data on a hardness-intensity diagram revealed that 
the whole Z-track was moving sideways forming a "paral- 
lel track" which is often seen. For this reason, we excluded 
these data from our analysis to avoid superimposing data 
in which the source has changed. Similarly, we excluded 
these data from the Z-track shown in Fig. 2. A spike in the 
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82.8/68 


5 


2.6±0.5 


1.99±0.16 


1.51±0.30 


2.4±0.5 


6.68±0.24 


6.15±0.45 


73.4/65 
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52.6/55 


Normal Branch 
















7 


2.6±0.5 


1.65±0.12 


1.43±0.50 


3.5±0.7 


5.95±0.16 


7.28±0.60 


78.4/68 
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7.9 ±1.3 
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0.9±0.9 


1.22±0.04 


4.80±1.30 


11.9±1.0 


5.53±0.60 


5.05±1.50 


47.9/47 
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for a distance of 10 kpc, the cut-off power law normalization is in units of photon cm" -2 s-i keV-i at 1 keV. 



data may also be seen at ^210 ksec lasting about 600 s, 
the nature of which is not clear, however, this event oc- 
curred in the data excluded from analysis. 

Spectra were extracted at 13 positions about equally 
spaced along the Z-track allowing for data in some parts of 
the horizontal branch being sparse. Narrow ranges of hard- 
ness ratio (0.015 wide) and intensity (30 count s"-'^ wide) 
were defined and incorporated in good time interval (GTI) 
files for each selection and used to extract PCA spectra 
and the corresponding background files. Deadtimc correc- 
tion was made using a local facility PCADEAD. Systematic 
errors of 1% were added to each channel; it was not nec- 
essary to regroup to a minimum count in each bin to al- 
low use of the statistic as the count in primitive bins 
was already high. A response function was generated using 

PCARSP. 

HEXTE spectra were extracted for each selection by 
using the same GTI files as for the PCA. The auxiliary 
response file (arf ) of May 2000 and the response matrix file 
(rmf) of March 1997 were used and the rmf was rebinned 
to match the actual number of channels in HEXTE spectra 
using the FTOOLS rddescr and rbnrmf. The maximum 
energies used in the PCA and HEXTE were set at the 
energies at which the source fiux became less than the 
background fiux, typically 30 keV in the PCA and 40 keV 
in HEXTE. 

The extended ADC emission model was applied in the 
spectral fitting in the form AB *(CPL + BB), where AB is 
absorption, CPL is a cut-off power law for Comptonization 
in an extended corona (Church & Balucinska-Church 
2004), and bb is the blackbody emission of the neutron 



star. It is well-known that the Comptonization cut-off en- 
ergy in the bright Z-track sources is relatively low at a few 
keV compared with higher values in less luminous sources. 
Thus there is a limited energy range for determination 
of the Comptonization power law index F and we have 
taken the approach of fixing T at 1.7, a physically rea- 
sonable value (Shapiro et al. 1976), as adopted previously 
for GX340-h0 and GX5-1 (Church et al. 2006; Jackson et 
al. 2009). When a best-fit was obtained, F was freed, but 
the value remained close to 1.7. In addition, a Fe line was 
detected at ~6.7 keV with an equivalent width of about 
70 eV which showed relatively little variation along the Z- 
track. Good fits were obtained to all spectra, and spectral 
fitting results are shown in Table 1 and in Figs. 3 and 4. 

In Fig. 3 we show results for the neutron star black- 
body temperature kT^B (upper panel) and the blackbody 
radius Rbb (lower panel) as functions of the total 1-30 
keV luminosity. A clear pattern of systematic variation is 
evident: the temperature is lowest at ^1.3 keV at the soft 
apex, suggesting that the mass accretion rate M is low. 
At this position Rbb is maximum at 12 km, consistent 
with the whole neutron star siirface being emitting. On 
the normal branch, the temperature rises while the black- 
body radius decreases, and the temperature continues to 
rise on the horizontal branch. 

Next in Fig. 4, we show the individual luminosities of 
the neutron star blackbody (Lbb) and the Comptonized 
emission (Xadc) as a function of the total 1-30 keV 
luminosity, so that neglecting the small contribution of 
the line: Lxot = -^bb + -^adc- It can be seen that the 
Comptonized emission is the dominant component, ten 



Balucinska-Church et al.: The nature of the Cyg-hke sources 



RXTE 



Cygnus X-2 



Horizontal Branch 




Soft Apex 



Flaring Branch ' Hard Apex 

. *' 

Normal Branch 



H 1 1 1 1 1 \ \ 1 1 1 1 1 H 



Flaring Branch 



Soft Apex 



Normal Branch 




Horizontal Branch 



Hard Apex 



15 



20 

Lp„t(10='ergs-') 



25 



Fig. 3. 



Left: blackbody temperature (upper panel) and radius 
(lower panel) as a function of the total luminosity. 
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Fig. 4. Luminosities of the continuum emission components: 
neutron star blackbody emisssion (upper panel) and ADC 
Comptonized emission (lower panel). 



times more luminous than the blackbody, and that as the 
source moves up the Z-track from the soft apex to the hard 
apex, this component doubles in luminosity. The X-ray 
intensity also, of course, increases by a similar factor and 
we suggest therefore that the mass accretion rate increases 
in this direction contrary to the widely-held view that M 
increases monotonically round the Z-track in the direction 
HB - NB - FB (Priedhorsky et al. 1986). These variations 
will be discussed in more detail in the next section in terms 
of all 3 Cygnus X-2 like sources, given their similarity of 
behaviour. 



4. The three sources: CygX-2, GX 340+0 and 
GX5-1 

We now present the results for Cygnus X-2 in compari- 
son with those we previously obtained for GX 340+0 and 
GX5-1 (Church et al. 2006; Jackson et al. 2009). Fig. 5 
(upper panel) shows the neutron star blackbody parame- 
ters. In all cases the temperature is lowest at the soft apex 
and rises continuously on the NB and HB. The blackbody 
radius at the soft apex in all sources lies between 10 and 12 
km (the dashed lines) suggesting that the whole neutron 
star is emitting. If this is assumed, the results provide a 
value of the neutron star radius based on the three mea- 
surements of 11.4 ± 0.6 km at 90% confidence. The flaring 



Cygnus X-2 



GX5-1 




L,„i (10^' erg s-') 



Fig. 5. Upper panel: blackbody temperature (upper panel) and 
radius (lower panel) as a function of the total luminosity; lower 
panel: the individual luminosities of the neutron star black- 
body (upper panel) and of the Comptonized ADC emission as 
a function of the total luminosity (lower panel). The soft apex 
(SA) and hard apex (HA) are marked. 



Balucinska-Church et al.: The nature of the Cyg-hke sources 



7 



branch is characterized by some increase of temperature 
and a marked increase of blackbody radius. 

The luminosities of the neutron star blackbody Lbb 
and of the Comptonized emission Ladc are shown in the 
lower panel. In the variation of Ladc with Lxot the three 
sources lie on a single line on the NB, suggesting that the 
sources differ only in their mass accretion rate, i.e. total 
luminosity. We make the standard assumption of accre- 
tion physics that the total luminosity of the source ixot 
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Fig. 6. Flux of the emitting part of the neutron star as a frac- 
tion of the Eddington flux (assumed constant): Z/Bdd/4 tt i?^, 
where R = 10 km. 



depends on the mass accretion rate M. A single line is, 
of course expected, since Labc comprises a large frac- 
tion of the total luminosity. The single line docs suggest 
that if, for example Cygnus X-2 became brighter it would 
appear similar to GX 340-1-0 in its X-ray properties. The 
increase of Labc in each source between the soft apex 
and the hard apex strongly suggests that the mass accre- 
tion rate M is increasing contrary to the standard view 
that M increases monotonically around the Z-track in the 
direction (HB - NB - FB). The observed increase of neu- 
tron star blackbody temperature on the normal branch is 
consistent with this as we would expect on the basis of 
simple accretion theory a higher temperature when more 
accretion flow reaches the surface of the neutron star. On 
the HB the luminosity of the Comptonized emission falls 
back towards its initial value. On the NB and HB there 
is an increase of blackbody temperature, but a decrease 
of blackbody radius. Consequently the blackbody luminos- 
ity increases less than expected for a temperature increase 
alone (~ T*). In flaring, there is an increase of fcTeB, al- 
though the flaring is weak in CygX-2. In all three sources, 
-^^ADC remains constant or falls in flaring, indicating that 
M does not increase. However, the blackbody luminosity 
increases substantially in the strong flaring in GX 340+0 
and GX5-1 without an increase in M. Thus the increase 
of Z/Tot is due to the increase in Lbb and this must clearly 



be caused by an additional energy source on the neutron 
star, i.e. unstable nuclear burning. 

4.1. Radiation pressure 

We next view the results obtained in terms of the radiation 
pressure of the neutron star which we can quantify based 
on the spectral fltting results. In all three sources fcTeB 
increases strongly on the NB and HB, by about a factor 
of two. Thus radiation pressure ^T"^ increases by nearly 
an order of magnitude. Moreover, as the source moves 
from the soft to the hard apex _Rbb decreases, so that 
the emission of the neutron star becomes concentrated, 
arising from an equatorial strip, from which the emissive 
flux becomes high. We consider the strength of the radi- 
ation pressure in terms of the ratio of the emissive flux / 
to the Eddington value /Edd (Sect. 1). Fig. 6 shows that 
this ratio is initially small at the soft apex (^-^0.2), i.e. sub- 
Eddington, but increases rapidly on the normal branch be- 
coming equal to unity at the hard apex, then continuing to 
increase on the horizontal branch having super-Eddington 
values as high as ///Edd = 3 in GX 340-1-0. The positions 
on the Z-track where the sources are super-Eddington are 
exactly those positions where radio emission is detected 
indicating the presence of jets, and so we propose that 
high radiation pressure plays a major role in launching 
the jets, by disrupting the inner disk and diverting ac- 
cretion flow into the vertical direction. In a more detailed 
discussion (Church et al. 2006), it is shown that the reduc- 
tion in blackbody radius on the NB and HB is consistent 
with this disruption. 

4.2. Unstable nuclear burning 

In Fig. 7, we show results derived from spectral fitting 
for the mass accretion rate per unit area of the neutron 
star TO, using M from the 1-30 keV source luminosity 
and dividing this by 47ri?33. The critical value tost at 
which there is a transition from stable nuclear burning 
(to > Tost) to unstable burning at smaller values of to 
(Bildsten 1998) is shown as a dashed line with its esti- 
mated 30% uncertainties (dotted). It may be seen that 
there is good agreement between the values of to at the 
soft apex and tost for the three sources suggesting that 
as a source descends from the hard apex towards the soft 
apex, miclear burning becomes ^instable and the result is 
the onset of flaring. On the flaring branch, the constancy 
of -Ladc (Fig- 5) indicates that the mass accretion rate is 
constant. The increase in source luminosity is due to an 
increase in Lbb at constant mass accretion rate and shows 
that there is an additional energy source on the neutron 
star not due to an increase of accretion rate, and nuclear 
burning is an obvious possibility. The above agreement 
with theory indicates that this is unstable nuclear burn- 
ing. This is further discussed in Sect. 6. During flaring 
there is an increase of blackbody radius (Fig. 5) from the 
value at the soft apex which we take to be the neutron star 
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radius. These increases, e.g. to 17 km in GX 340+0, indi- 
cate an expansion of the unstable nuclear burning region 
beyond the surface of the neutron star similar to radius 
expansion in the burst sources. 
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Fig. 7. Mass accretion rate per unit neutron star emitting area 
rh. 



5. The kHz QPO 

We next present results on the kHz QPO in Cygnus X-2 
and GX 340+0 and compare these with our previous re- 
sults for GX5-1 (Jackson et al. 2009). The approach 
adopted for GX5-1 was to use the same selections of 
data made around the Z-track for both spectral analy- 
sis and timing analysis, thus allowing the QPO results to 
be viewed in terms of the spectral and physical evolution 
on the Z-track. In the case of Cygnus X-2, the same gen- 
eral approach was used, however, because of the relative 
weakness of the kHz QPO, it was necessary to increase 
the size of the selection boxes (from those in Fig. 2) to 
allow detection of the QPO. Spectral fitting was again 
carried out for the new, larger boxes. The results of this 
were, however, entirely consistent with the results of Sect. 
3, and were used in analysis of the QPO results. In Fig. 
8 we show the new selection boxes labelled from 1 to 10. 
As is well-known, kHz QPO are only detected on the HB, 
at the hard apex and sometimes at the upper part of the 
NB, so the selection numbering begins at the end of the 
Z-track where kHz QPO are expected. 

In the case of GX 340+0, timing analysis of the ob- 
servation of 1997 September for which we previously pre- 
sented spectral fitting results (Church et al. 2006) revealed 
the presence of kHz QPO on the HB. However, because of 
a sparsity of data in the centre of the horizontal branch, 
a sufficiently good dependence of kHz QPO frequency on 
branch position could not be obtained. Consequently the 
observation of 1997 June 6-10 was also analysed which 
was reasonably close in time to the original data. Spectral 
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Fig. 8. Selection boxes used for fitting of the power spectra. 
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Fig. 9. Kilohortz QPO frequencies as a function of source lu- 
minosity for Cygnus X-2, GX 340+0 and GX5-1. 



analysis as a function of Z-track position produced spec- 
tral fitting results very similar to those presented in Figs. 
5, 6 and 7. Better coverage of the HB was available, and 
we used this observation carrying out both spectral and 
timing analysis as a function of position on the Z-track. 

In both CygX-2 and GX 340+0, power spectra were 
extracted as described in Sect. 2 and fitted in the range 150 
- 4096 Hz, using a model consisting of two Lorentzian lines 
for the kHz QPO plus a power law which was able to fit the 
continuum adequately within this frequency range. KHz 
QPO could be detected on the horizontal branch only; in 
Cyg X-2 only the higher frequency kHz QPO could be de- 
tected in box 1 (at the end of the Z-track) . Figure 9 shows 
the best-fit lower QPO [vx) and upper QPO (1/2) frequen- 
cies as a function of the total 1-30 keV luminosity which 
was obtained by spectral fitting using the same model as 
in Sect. 4. The results for Cyg X-2 and GX 340+0 are com- 
pared with the corresponding results for GX5-1 (Jackson 
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et al. 2009) also shown in Fig. 9. All three sources dis- 
play an increase of QPO frequency with luminosity during 
movement along the horizontal branch as is well-known 
from previous work (e.g. Jonker et al. 2000; Wijnands et 
al. 1998a, b). The main aim is, of course, to find the cause 
of the frequency change. 

The frequency V2 is known to correspond to the or- 
bital frequency in the accretion disk at positions a little 
way into the disk and we next view the results in these 
terms. In each case we derive the radial position r at which 
the Keplerian orbital frequency equals the measured QPO 
frequency v, using the relation 



277 1/ : 



GM 



Thus the radius r decreases with increasing Lxot for both 
upper and lower QPO. We show in Fig. 10 the radial po- 
sition r not as a function of luminosity, but as a function 
of the ratio ///Edd which is a measure of the strength 
of the radiation pressure close to the neutron star (Sect. 
4.1). The figure includes our results previously published 
for GX5-1 (Jackson et al. 2009). In the case of GX5-1, 
a strong increase of the QPO radius took place when the 
increasing blackbody temperature moving from the hard 
apex on the HB caused a marked increase in ///Edd at 
about / = /Edd- Super-Eddington values of the emissive 
flux suggest disruption of the inner disk which will in- 
crease with ///Edd and if the upper frequency oscillation 
takes place at the inner disk edge, the radius r will move 
outwards as observed. 
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Fig. 10. Keplerian radial positions corresponding to the upper 
kHz QPO {v2) and lower kHz QPO {v\) as a function of the 
strength of the neutron star radiation pressure expressed as 
///Edd (see text). 



In CygnusX-2 and GX 340-1-0, the radial position r 

also increases with ///Edd as in GX5-1. However, in these 
sources kHz QPO could only be detected on the horizontal 
branch so there are no QPO detections for ///Edd < 1- 



In these sources, it appears that higher values of ///Edd 
are needed to produce a substantial change in QPO ra- 
dius, and a possible reason for this is discussed in Sect. 
6. However, the results do suggest an explanation of the 
kHz QPO and their variation along the HB based on our 
spectral fitting results. In this explanation the disk is dis- 
rupted on the upper normal and horizontal branches as 
the neutron star emissive flux becomes super-Eddington. 
The higher frequency QPO with frequency (having the 
smaller radial position) is at all times an oscillation on 
the inner edge of the disk which moves outwards with in- 
creasing radiation pressure and increasing disk disruption, 
causing the change in the measured QPO frequency. The 
lower frequency QPO {vi) corresponding at the hard apex 
to r ~20 km also moves to much larger radial positions 
~30 - 40 km (Fig. 10). The mechanism producing this os- 
cillation is not clear, but as its variation follows that of 
1/2, it seems clear that its frequency is derived from, i.e. 
depends on, V2 

6. Discussion 

We have shown that application of the extended ADC 
emission model for LMXB provides good fits to the spectra 
of all three Cygnus X-2 like Z-track sources at all positions 
along the Z-track. Moreover, the physical interpretation 
of the results is straightforward and suggests a plausible 
explanation of the Z-track phenomenon in this group of 
Z-track sources. Our explanation of the Z-track is that 
the soft apex is the lowest luminosity state of the source 
with a low value of M, with emission taking place from 
the whole neutron star which has its lowest temperature. 
On the normal branch, the increase of intensity and ADC 
luminosity suggest an increase of M leading to heating of 
the neutron star and a strong increase in radiation pres- 
sure close to the neutron star as shown by the increase 
of ///Edd (Fig. 6), which continue on the HB. We sug- 
gest that this has a strong effect on the inner accretion 
disk causing disruption of the disk. The horizontal force 
may not directly remove matter from the disk, but be- 
cause the unperturbed height of the inner disk in LMXB 
at these luminosities greater than 10^* erg s~^ is 20 - 50 
km, the radiation pressure can also act in a direction close 
to vertical blowing away material from the upper layers of 
the disk. For the strongly super-Eddington fluxes that we 
measure close to the equatorial emitting zone of the neu- 
tron star, a fraction of the mass accretion flow may be 
diverted vertically upwards and be ejected from the sys- 
tem as massive blobs of plasma forming the jets above 
and below the disk. Thus we propose that high radia- 
tion pressure is a necessary condition for jet formation. 
There are large differences in luminosity between the Z- 
track sources with Cygnus X-2 being typically 4 times less 
luminous than ScoX-1. Although the sources are detected 
in radio, in Cygnus X-2 jets have not so far been imaged 
as in ScoX-1 (Fomalont et al. 2001), so there is the possi- 
bility that a second requirement for strong jet formation 
is a higher value of M. 
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Table 2. The four regimes of stable or unstable nuclear burn- 
ing from Bildstcn (1998), showing the critical values of rfi that 

demarcate the regimes. 



12 3 4 

CNO Hot CNO He unstable H,He 
unstable stable in H-rich stable 

environment 

Critical values of m: 
1 X 10''' g cm^'^ s^^ (regimes 1-2) 
5 X 10^ g cm"^ s"^ (regimes 2 - 3) 
1.3 X 10^ g cm"^ s"^ (regimes 3 - 4) 



The results demonstrate that in all of the three 
CygnusX-2 like Z-track sources, the flaring branch con- 
sists of a strong increase of neutron star blackbody lumi- 
nosity while the Comptonized emission luminosity is con- 
stant (or decreases slightly) showing that M is constant. 
This indicates that the FB consists of unstable nuclear 
burning on the neutron star. In the theory of unstable 
nuclear burning (Fujimoto et al. 1981; Fushiki & Lamb 
1987; Bildsten 1998; Schatz et al. 1999) the physical con- 
ditions in the atmosphere of the neutron star depend on 
the mass accretion rate per unit area m, i.e. M divided 
by the emitting area. The various regimes of stable and 
unstable burning are shown in Table 2. The theory applies 
to a wide range of mass accretion rates, although it was 
developed to describe X-ray bursting at lower mass accre- 
tion rates, and also applies to high mass accretion rates 
typical of the Z-track sources, and the transition between 
regimes 3 and 4 takes place for luminosities of about the 
Eddington limit. In regime 3, burning is unstable He burn- 
ing in a mixed H/He environment (Bildsten 1998) but for 
rh > rhsT, H and He burning of the accretion flow become 
stable. The agreement in all three sources of the mea- 
sured rh values (Fig. 7) with the critical value within its 
uncertainties clearly supports flaring being unstable nu- 
clear burning. The value of to at the soft apex increases 
systematically with Lxot suggesting that other, secondary 
effects may take place. However, the results thus provide 
evidence in two ways that the FB in all three CygnusX-2 
like Z-track sources consists of unstable nuclear burning. 
There is also an expansion of the emitting region beyond 
the surface of the neutron star and we propose that radius 
expansion may take place as it does in some X-ray bursts 
in lower luminosity sources. 

The proposed explanation of the Z-track in the 
Cygnus X-2 like sources involves an increase of M between 
the soft apex and hard apex, and a constant M in flaring. 
Thus we note that this explanation departs from the often 
held view that the mass accretion rate increases monoton- 
ically around the Z-track in the direction HB - NB - FB. 

Wo arc also carrying out analysis of ScoX-1 like 
sources aimed at discovering the fundamental differences 
between the Cyg-like and Sco-like groups. Preliminary 



work on the source GX 17-1-2 shows that the behaviour 
on the normal and horizontal branches is similar to that 
of the Cygnus X-2 like sources. However, in flaring there is 
clear evidence for an increase of mass accretion rate. Thus 
it appears that on both the FB and NB, M increases as we 
move away from the soft apex, however with rather differ- 
ent conditions on the neutron star, i.e. different values of 
kT and Rbb- Clearly, to fully understand the ScoX-1 like 
sources extensive work will be required; however, it does 
appear that the difference between the two groups relates 
to the flaring branch, and not, for example, to inclination. 

Relevant to this is the decay of the transient LMXB 
XTE J1701-462, which foUowing an outburst in 2006 
January, decreased by more than an order of magnitude in 
luminosity over about 500 days. Lin et al. (2009a) showed 
that the source made a transition from a Z-track source 
of ScoX-1 like type to Atoll behaviour, these changes be- 
ing apparent from hardness-intensity diagrams. Towards 
the end of the decay, three X-ray bursts were observed 
and Lin et al. (2009b) claim that their results cast doubt 
on our association of the FB with unstable nuclear burn- 
ing. However, the claim is not valid. Two of the bursts 
occurred in the time period labelled "V" by Lin et al. in 
which the source showed clear Atoll behaviour, having a 
banana branch in hardness-intensity. The other burst how- 
ever occurred in period "IV" just before period V began, 
when flaring was weak and which in period V disappeared 
altogether. Although the luminosity was at that time 10% 
of Eddington, the authors describe the source at the pe- 
riod of this burst as being on the FB of a Z-track state. 
On the basis that it may be unlikely that bursting and 
unstable nuclear burning in flaring take place at the same 
time it was claimed that our association of the FB with 
unstable burning must be doubted. The argument is not 
correct (Lin et al. 2009a, b) as our results and model for the 
Z-track specifically refer to the Cygnus X-2 like sources. 
Thus the observation of the burst during a period when 
the source was described as being in a ScoX-1 like state 
(but normally described at such low luminosity as an up- 
per banana state) does not cast doubt on our explanation 
of the Cygnus X-2 like sources in which the evidence that 
the FB consists of unstable burning is strong. 

The kHz QPO results support the idea that high ra- 
diation pressure of the neutron star on the upper NB and 
HB disrupts the inner disk and that the upper QPO at 
frequency 1/2 is an oscillation actually at the inner edge 
of the disk. The inner radius is variable because of in- 
creasing disruption on the horizontal branch thus causing 
the well-known variation of 1^2 on that branch. The fact 
that the radial position associated with the QPO increases 
when there is also an increase of ///Edd is thus sugges- 
tive. Close to the soft apex, ///Edd is much less than unity 
and we do not expect the disk to be disrupted at all. Thus 
the inner edge of the disk will touch the neutron star, so 
that oscillations of the edge will probably not be possi- 
ble. Thus this model for the upper frequency kHz QPO 
predicts that kHz QPO will not be detected on the lower 
normal branch, in agreement with observation. 
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Fig. 11. Inner disk height profiles H{r) for the three sources 
with points indicating the position of the inner disk edge in 
GX5-1 and CygnusX-2 for ///Edd = 1-3 (see text). 



In GX5-1, the most luminous source, substantial 
change in the radial position of both upper and lower 
QPO takes place at ///Edd ^1, but in GX 340+0 and 
CygnusX-2, larger values of ///Edd are needed which may 
appear surprising. However, the inner disks of the three 
sources have very different vertical heights and height pro- 
files H{r) and these are shown in Fig. 11 based on the 
standard theory of the thick disk supported vertically by 
its own radiation pressure (Frank et al. 2002). For each 
source, H{r) is plotted to correspond to the luminosity 
at the hard apex of the Z-track where we expect there 
to be disruption of the inner disk. At radial positions r 
< 50 km, appropriate to the kHz QPO, H{r) is rising 
and eventually reaches the equilibrium height of the ra- 
diatively supported disk. Also marked on the figure arc 
radial positions corresponding to a fixed value of ///Edd 
= 1.3. Fig. 10 shows that for this value, the upper kHz 
QPO has an equivalent orbital radius of 24 km GX5-1 
and 17 km in CygX-2 (there is no value of ///sdd that 
includes all three sources). Thus the two points plotted 
in Fig. 10 indicate the radial positions where the QPOs 
occur, which we assume are the inner edges of the disk 
and so indicate the radial and vertical extent of disk dis- 
ruption. This shows that in the brighter source GX5-1 
the disruption for a given level of radiation pressure is 
greater. A theoretical treatment of disruption is complex 
and beyond the scope of this work; however, it does not 
appear unreasonable that in a source with greater disk 
height (GX5-1), the disruption is greater. 
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